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Abstract 

We consider the incompressible Navier-Stokes (NS) equations on a torus, 
in the setting of the spaces 1? and H 1 ; our approach is based on a general 
framework for semi- or quasi-linear parabolic equations proposed in the pre- 
vious work [9]. We present some estimates on the linear semigroup generated 
by the Laplacian and on the quadratic NS nonlinearity; these are fully quan- 
titative, i.e., all the constants appearing therein are given explicitly. 
As an application we show that, on a three dimensional torus T 3 , the (mild) 
solution of the NS Cauchy problem is global for each H 1 initial datum uq with 
zero mean, such that ||curluo||^2 ^ 0.407; this improves the bound for global 
existence ||curluo||.L 2 ^ 0.00724, derived recently by Robinson and Sadowski 
[10]. We announce some future applications, based again on the H 1 framework 
and on the general scheme of [9]. 
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1 Introduction 



As well known, the great open problem about the incompressible Navier-Stokes (NS) 
equations is to prove global existence of the solutions in three space dimensions, for 
all sufficiently smooth initial data. 

A much more modest, but realistic research program about these equations is the 
following. 

(i) Deriving estimates for the linear semigroup generated by the Laplacian and for 
the NS quadratic nonlinearity, within an appropriate functional setting (say, the 
Sobolev spaces on R 3 or on the torus T 3 ). This should be done paying attention to 
the strictly quantitative aspects, and having as a goal the best achievable accuracy. 

(ii) Explicitating the estimates for the NS solutions derivable from the previous 
framework. 

Some typical issues to be treated as parts of item (ii), always using the information 
from (i), are the following ones: 

(ii.a) Deriving quantitative lower bounds on the interval of existence of the solution, 
for a given initial datum. 

(ii.b) Producing a sufficient condition for global existence of the solution, when the 
norm of the initial datum is below an explicitly given upper bound, 
(ii.c) Devising a posteriori tests on the approximate solutions (e.g., of Galerkin type) 
to get quantitative estimates on exact solutions. 

A number of papers published in recent times are more or less related to the above 
research program. We mention, in particular: the works by Chernyshenko, Con- 
stantin, Robinson, Titi [1], Dashti and Robinson [2] on the approximate solutions 
of the NS equations, and on the a posteriori estimates derivable from them; the 
paper by Robinson and Sadowski [10] on the conditions for the existence of (exact) 
global solutions, containing very interesting considerations on the computational 
times required to check them via suitable approximate solutions; our papers [7] [9] 
on approximate solutions and a posteriori estimates for general semi- or quasi-linear 
evolution equations, including (in the second reference) some applications to NS 
equations on a torus. Let us also mention that a similar attitude towards approx- 
imate solutions and a posteriori bounds has been developed by Machiels, Peraire 
and Patera [6], Veroy and Patera [12] in connection with the NS equations (or some 
related space discretizations), dealing mainly with two-dimensional steady problems. 

The present work outlines a framework to treat problems (ii.a) (ii.b) (ii.c) for the 
incompressible NS equations on the torus T 3 . The main application presented here 
is a fully quantitative bound on the initial datum yielding global existence, in the 
spirit of (ii.b); a few words will be spent on other issues within the scheme of (ii), 
to be treated elsewhere. 

Our present constructions are related to the setup of [9], where the basic functional 
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space for the three dimensional NS equations was the Sobolev space H" (T 3 ) C 
H"(T 3 ), with n > 3/2 (the subscripts E, indicating the conditions of zero diver- 
gence and zero mean, respectively). 

On the contrary, here the solutions of the NS equations take values in H^ () (T 3 ) C 
Hg(T 3 ). EIq(T 3 ) is equipped with the norm v h-> \\v\\i := \\V— A v\\ L 2 (A the Lapla- 
cian); the latter is restricted to H^ (T 3 ), where it is found to coincide with the L 2 
norm of the vorticity: 

\\v\\i = 1 1 curl i; 1 1 x,a . (1.1) 

Throughout the paper, the incompressible NS equations for the velocity field u are 
written as 

(in 

^f = Au- Z{u.du) , (1.2) 
at 

with £ indicating the Leray projection on the divergence free vector fields. As in 
our previous work, we emphasize the role played by the constants in certain basic 
estimates on the semigroup (e iA ) and on the NS bilinear map (v,w) i— > £(v»dw). In 
particular, for 1/2 < uo < 1 we consider the negative order Sobolev space EI~^(T 3 ) 
with the norm || • ||_ w := ||\/A ■ ||^2, and give explicitly a constant K w such that, 
for all v,w G H^T 3 ), 

\\£,(v.dw)\\- u ^/r w |H|i||iw||i ; (1.3) 

for example, if uo = 7/10 the above inequality is fulfilled with K 7 / 10 = 0.361. On 
the grounds of our previous investigations on multiplication in Sobolev spaces [8], 
the constant K u determined for any uo by our method is probably close to the best 
constant fulfilling (1.3). 

The evaluation of K u can be combined with some strictly quantitative estimates on 
(e tA ), and with the general setting of [7], to obtain a number of results on the exact 
NS solutions. In particular, here we show that the incompressible NS equations have 
a global solution u for any initial datum u G H^ (T 3 ) such that 

||«o||i < 0.407; (1.4) 

this improves the quantitative condition for global existence presented in [10] (page 
43) that, in our notations, can be written as ||wo||i ^ 0.00724. In connection with 
these results, let us mention that the idea of an explicit norm bound, sufficient 
for global existence, dates back to the seminal paper by Kato and Fujita [3]; here, 
an approach slightly different from ours was outlined but not developed in a fully 
quantitative way, since the authors did not compute all the necessary constants (nor 
indicated how to do the missing calculations). 

Let us outline the organization of the present paper. In Section 2, we review some 
part of the abstract framework of [9] for semi- or quasi-linear evolution equations; we 
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consider in particular the case of a quadratic nonlinearity, suitable for application 
to the NS equations. 

In Section 3, we fix our standards about Sobolev spaces on a torus T d (for the mo- 
ment, of any dimension d ^ 2), and indicate the connections of this formalism with 
the NS equations (of course a precise definition of these standards, and especially 
of the Sobolev norms, is required to make meaningful the quantitative estimates 
mentioned before, to have comparisons with [10] and other papers, etc.). 
In Section 4 we discuss the semigroup (t, v) \— > e tA v as a map [0, +oo) x H^ (T d ) — > 
H^ (T d ), or as a map (0, +oo) x H.~"(T d ) -> H^ () (T d ), for d ^ 2 and < u < 1. 
In Section 5, after some preliminaries we discuss the bilinear map HI* (T d ) x H^ (T d ) 
— > H~ w (T <i ), (v,w) i— > £(v»dw) for d > 2 and a; > d/2 — 1; we infer an inequality 
of the type (1.3) and show how to compute the related constant K w (in fact, also 
dependent on the dimension d). 

In Section 6, we put together the results of Sections 4-5 on (e tA ) and on the NS 
bilinear map, which is possible if (d = 2, < u < 1 or) d = 3, 1/2 < uj < 1; this 
gives our final "if 1 framework" for the NS equations on T 3 . As an output of this 
framework, choosing u = 7/10 we obtain the condition (1.4) of global existence in 
H^ (T 3 ), together with a quantitative estimate on the exponential decay rate of the 
global solution u. 

In Section 7 we present a method to get estimates on the (exact) NS solution ana- 
lyzing a posteriori any approximate solution, in the H 1 framework on T 3 . Here we 
specialize the general method of the "control inequality" introduced in [7] [9]; this 
was already applied in [9] to the NS equations on H™ (T 3 ), n > 3/2, for example 
to get estimates from the Galerkin approximate solutions. Some applications of the 
control inequality to the approximate NS solutions in H^ (T 3 ) will be presented 
elsewhere. 

The paper is completed by some Appendices, concerning more technical issues. In 
Appendix A, we prove some estimates on the semigroup (e* A ) employed in Sections 
3-4. In Appendix B we illustrate some general facts on (discrete) convolutions 
of unimodal functions; these results are employed in the subsequent Appendix C, 
where we study a convolution whose maximum gives the constant K u of Eq. (1.3). 
In Appendix D we give some details related to the computation of K u in dimension 
d = 3, for u = 7/10. 

2 An abstract framework for evolution equations 
with a quadratic nonlinearity 

In [9], we considered a framework for evolutionary problems (mainly, of parabolic 
type) with a nonlinear part of a fairly general kind; in the cited paper, the framework 
was subsequently specialized to the case of a quadratic nonlinearity. For our present 
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purposes, it suffices to review the case of a quadratic, time independent nonlinearity. 
The framework. Throughout this section, we consider a set 

(F + ,F,F_,„4,J>) (2.1) 

with the following properties. 

(PI) F + , F and F are Banach spaces with norms || || + , || || and || ||_, such that 

F + ^ F ^ F (2.2) 

(the symbol <^-> indicating that one space is a dense linear subspace of the other, and 
that the natural inclusion is continuous). Elements of these spaces are generically 
denoted with v, w, ... . 
(P2) A is a linear operator such that 

A : F + -> F_ , v i Av . (2.3) 

Viewing F + as a subspace of F_, the norm || || + is equivalent to the graph norm 
i> G F + i— > ||d||_ + 

(P3) Viewing A as a densely defined linear operator in F_, it is assumed that A 
generates a strongly continuous semigroup (e tA ) te \o j+00 ) on F (of course, from the 
standard theory of linear semigroups, we have e tA (F + ) C F + for all t ^ 0). 

(P4) One has 

e tA (F) C F for t G [0, +oo) ; (2.4) 

the function (v,t) i— > e tA v is continuous from F x [0, +oo) to F, yielding a strongly 
continuous semigroup on F as well. Furthermore, there is a constant B > such 
that 

lle'^H < e~ B *|H| forall^GF. (2.5) 

(P5) One has 

e tA (F_) C F for t G (0, +oo) ; (2.6) 

the function (v,t) h- > e tA v is continuous from F_ x (0, +oo) to F (in a few words: 
for all t > 0, e tA regularizes the vectors of F_, sending them into F continuously). 
Furthermore, there is a function \i G C((0, +oo), (0, +oo)) such that 

lle'^H ^ fi(t)e' Bt \\v\\_ for t > 0, v G F_, B as in (2.5) . (2.7) 

The function /x behaves like an integrahle power of 1/t for t — > 0, i.e., 

/i(t) = 0(l/t e ) for i -> 0+ , £ G [0, 1) . (2.8) 

Finally, there is a constant N G (0, +oo) such that 

/ ds e~ Bs fi(t - s) ^ N for t G [0, +oo) . (2.9) 
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(P6) CP is a bilinear map such that 

CP : F x F -> F_ , (u, w) h-> 0>(u, w) ; (2.10) 

we assume continuity of CP, which is equivalent to the existence of a constant K G 
[0, +oo) such that, for all v, w G F, 

||CP(v,w)||_ ^ #WHI • (2.11) 

The initial value problem; some results of uniqueness and existence. The 

initial value problem with initial datum uq G F is the following one: 

Find u G C([0, T),F) such that 

u(t) = e tA u + f ds e {t - s)Af ?(u(s),u(s)) for all t G [0, T) . (2.12) 
Jo 

Of course, we say that a solution u is maximal if it cannot be extended and global 
if its domain is [0, +oo). The following facts are well known. 

(i) With a bit more regularity (u G F + and F reflexive), (2.12) is equivalent to 
the Cauchy problem du/dt = Au + CP(ii, u), u(0) = Uo for an unknown function 
u G C([0, T),F + ) fl C 1 ([0, T), F_). Independently of these stronger assumptions, a 
solution u of (2.12) is usually called a "mild solution" of the Cauchy problem. 

(ii) (2.12) has a unique maximal solution (and any other solution is a restriction of 
the maximal one); 

(iii) (2.12) has a global solution if the initial datum is small. Among the many ref- 
erences available on this point, for convenience we refer to our work [9], Proposition 
5.12, yielding the following statement: 

2.1 Proposition. Suppose 

AKN\\u \\ ^ 1 ; (2.13) 
then, the problem (2.12) has a global solution u : [0, +oo) — > F. Furthermore, 

\\u(t)\\ < X(4KN\\u \\) e- Bt \\u Q \\ for t G [0, +oo) , (2.14) 
where X G C([0, 1], [1,2]) is the increasing function defined by 

X(z) := I 1 ~$q~ Z A"" ^(0,1], (215) 
[ 1 for z = . 
More roughly, due to the features of the function X , one has 

\\u(t)\\ < 2e~ Bt \\u \\ for t G [0, +oo) . (2.16) 

Let us mention that the global solution u in Proposition 2.1 is obtained as the 
limit of a Picard iteration: u(t) = limfc_» +00 Uk(t) for each t G [0, +oo), where 
(uk)k=o,i,2,... is the sequence of functions in C([0, +oo), F)) defined by u (t) := and 
u k+1 (t):=e tA u + f*dse^ A ?(u k (s),u k {s)). 
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3 Sobolev spaces and the Navier-Stokes equations 
on a torus 



In this section we consider any space dimension 

d^2. (3.1) 

We use r, s as indices running from 1 to d; elements a, b,.. of H d or C d are written 
with upper or lower indices, according to convenience: (a r ) or (a r ), (b r ) or (b r ). For 
a, b G C d (say, with upper indices), we put 

d 

a»b := S ^a r b r \ \a\ := \Ja»a (3.2) 

r=l 

where a := (a r ) is the complex conjugate of a. Hereafter we refer to the d- 
dimensional torus 

T := R/(2vrZ) , (3.3) 




whose elements are typically written x = (^ r ) r =i,...d- 

Distributions on T d , Fourier series and Sobolev spaces. We introduce the 
space of periodic distributions D'(T d , C) = D' c , which is the dual of C°°(T d , C) = 
(equipping the latter with the topology of uniform convergence of all deriva- 
tives); (v, f) denotes the action of a distribution v G D' c on a test function / G C^. 
We also consider the lattice Z d of elements k = (/c r ) r= i i ... i d and the Fourier basis 
(ek)kez", where 

e*:T"-C, ek (z):=j^e*-* (3.4) 

(k»x = J2t=i^rX r makes sense as an element of T d ). Each v G D' c has a unique 
(weakly convergent) Fourier series expansion 

v Vk6k ' Vk := ( u ' e -fc) G C • ( 3 - 5 ) 

kez d 

(As well known, the Fourier coefficients v k of any v G D' c grow polynomially in \k\\ 
conversely, any family (vk)k&z d of complex numbers with such a polynomial growth 
is the family of the Fourier coefficients of some v G D' c ). The mean of v G D' c is 
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(of course, (v, 1) = f Td dxv(x) if v G L 1 (T d )); the space of zero mean distributions 



is 



D' co ■= {v G D' c I (v) = 0} . (3.7) 

The distributions v in this subspace are characterized by the equivalent condition 
t>o = 0; so, their relevant Fourier coefficients are labeled by the set 

Z^:=Z d \{0}. (3.8) 

The complex conjugate of a distribution v G D' c is the unique distribution v such 
that (v, f) = (v, f) for each /eC"; one has v = Y^kez^^k- 
The distributional derivatives d/dx s = d s (s — l,...,d) and the Laplacian A := 
z\2t=i @ss obviously send D' c in D' co , and are such that, for any v G D c , d s v = 
i ^2 k€Z d k s Vk&k, At> = — X^fcezg \k\ 2v k&k- For any n G R, we further define 

V=A n : D' c - D' co , u ^ v^A\ := ^ |£f ^e, . (3.9) 

fcez^ 

In the sequel, we are interested in the space of real distributions D'(T d , R) = D', 
defined as follows: 

D' := {v G D' c | v = v} = {v G D' c \ = v^ k for all k G Z d } ; (3.10) 

of course, v E D' implies (v) G R. We also set 

^:={^5' (u) = 0}| (3.11) 

all the differential operators mentioned before send D' into D' . 

Hereafter we consider the real Hilbert space L 2 (T d ,R,dx) = L 2 , with the inner 

product (v \w)l 2 '■= f Td v(x)w(x)dx = J2 k&Z d ^k w k an d the induced norm 



MU 2 = \ / ^ 2 (a:)^ = /J^ |w fc | 2 . 
V jTd V fcez d 



(3.12) 



To go on, we introduce the zero mean Soholev spaces H™(T d , R) = H™. For each 
n G R, 

H£ := {v G D ( ' | ^f^K n v G L 2 } = {«£ D' \ l*f "M" < +oo } ; (3.13) 

this is a real Hilbert space with inner product (v\w) n := (\ / —A n y | V-A"^)^ 
= X^fcez d l^| 2n W^fc and the induced norm 



\v\\ n =\\y/=A n v\\ L 2= /]T|A;| 2 «M 2 . (3.14) 



Let us consider, in particular, the case n is a nonnegative integer. For all v G 
D> and k G Z d , one has \k\ 2n \v k \ 2 = (£j=i*2)>fc| 2 = E£,...,. B =i ^-^J^l 2 = 
El,..., Sn =i l(9 ai ... s „u) fc | 2 ; so, 

# " := {v G Z^, I 9 s ,.. Sn ^ G L 2 Vs 1? s n G {1, d}} ; (3.15) 



A W d ^-Sn V \ 
\ Sl,...,S„ = l 



2 

L 2 



Vu G # n . 



(3.16) 



All these details about Sobolev spaces may seem pleonastic, but in fact they are 
useful to ensure that some norms we use later coincide exactly with the norms of 
[10]. 

For (integer or noninteger) 



n 



< n, one has HI 



and 



U on H£ 



In particular, H° is the subspace of L 2 made of zero mean elements. For any real 
n, Ai7™ = H'q~ 2 and A is continuous between these spaces. Obviously enough, we 
could define as well the complex Hilbert spaces L 2 and H™ ; however, these are 
never needed in the sequel. Other facts about Sobolev spaces (in particular, the 
duality between H~ n and if") are mentioned when necessary in the sequel. 

Spaces of vector valued functions on T d . If V(T d , R) = V is any vector space 
of real functions or distributions on T d , we write 



Y{T d ) = V := {v = (v\ ...,v d ) \v r eV for all r} . 



(3.17) 



In this way we can define, e.g., the spaces W(T d ) = W, h 2 (T d ) = L 2 , H"(T d ) = H™. 
Any v = (v r ) G D' is referred to as a (distributional) vector Geld on T d . We note 
that v has a unique Fourier series expansion (3.5) with coefficients 



v k = K)r=i,...,<i e C c 



(3.18) 



as in the scalar case, the reality of v ensures = v_ k . We define componentwisely 
the mean (v) G R^ of any v G D' (see Eq. (3.6)); D' is the space of zero mean 
vector fields. We similarly define componentwisely the operators d s , A, \j — A : 
D' — > D',. L 2 is a real Hilbert space; its inner product is as in the line before 
(3.12), with v(x)w(x) and VkW k replaced by v(x)»w(x) = Er=i yr { x ) w r (x) and 
Uk» Wk = E r =i v k w l- For an y rea ^ n -> ^ ne n-th Sobolev space of zero mean vector fields 
H"(T d ) = H™ is made of all d-uples v with components v r G if"; an equivalent defi- 
nition can be given via Eq.(3.13), replacing therein L 2 with L 2 . H" is a real Hilbert 



space with the inner product (v\w) n := (\ / —A n v \ \ / —A n w) L 2 = J2kezf t 1^ 



2n 



the induced norm || || n is given, verbatim, by Eq. (3.14). If w G H™ has components 
w r , we have obviously 



W r, = 



n ^ 



E 

r=l 



r\\2 . 

n i 



(3.19) 
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for n a nonnegative integer, the vector analogue of Eq. (3.16) is 



\w\ 



v E W^-s^ r \\h ■ (3-20) 

\ r,si, 



Divergence free vector fields. Let us consider the divergence operator div : 
D' — > Dq, f I— > divt> := Ylt=i 9 r v r ] of course, divt> = i ^2 keZ ^* v k) e k- The space 
of divergence free (or solenoidal) vector fields is 

D' E := {v G D' | divv = 0} = {v G O' | fc.v fe = VA; G Z d } ; (3.21) 

in the sequel, we consider as well the spaces 

O' E0 :=D' E nD;, IL:=»' E nH; (n G R) (3.22) 

(as usually, reference to T d is omitted for brevity: for example, H E0 stands for 
H E0 (T d )). H E0 is a closed subspace of the Hilbert space H™, equipped with the 
restriction of the inner product ( | ) n . We note the relations A(B' E ) = D' E0 and 
A(H E +2 ) = H E0 , following immediately from the Fourier representations. 
The Leray projection is the map 

£ : W -> D' E , v ^ £v := ^ (£ fc ^H , (3.23) 

feez d 

where, for each fc, £fc is the orthogonal projection of C d onto the orthogonal com- 
plement of k] more explicitly, if c G C d , 

£ c = c , £ fc c = c - ^ k for fceZf . (3.24) 

From the Fourier representations of £, ( ), <9 S , etc., one easily infers the following 
statements: 

(£u) = (u), £(<9 s w) = d s (£u), fi(Au) = A(£u) Vu G D'; (3.25) 

£D' = D E0 , £H" = e E0 ; ||£u|| n ^ ||u|| n Vn G R, v G H£ . (3.26) 

The spaces H*, H| . In the sequel, we are often interested in specializing the 
previous considerations to the case n — 1. This carries us to the space 



Hj := G D' | v 73 ^ v G L 2 } = {ve O' | <9X G L 2 Vr, s G {1, d}} 

= {^eD;,| ^ |A;| 2 M 2 < +oo } , (3.27) 



with the norm 

\\v ||i = || V— A v\\ L 2 = 



\ r,s=l 



We are as well interested in the divergence free part H B0 = D' E n Hj. We note that 

||u||i = ||curlv|| L 2 if d = 3 and v G H E0 , (3.29) 

where the right hand side contains the curl of vector fields (i.e., (curlf) 1 := d 2 v 3 — 
83V 2 , etc., intending derivatives in the distributional sense). To check Eq. (3.29), one 
notes that (curlv)^ = %k A v\. where A is the (complexified) vector product (so that 
(k A a) 1 := k 2 a 3 — k^a 2 , etc., for all a G C 3 ); if v is divergence free, the orthogonality 
property k»v k = implies \k A Vk\ = \k\\vk\, whence ||curlv||^ 2 = J^feez 3 I^Pl^fcl 2 = 
\Ml 

(Before going on, we point out that the norm ||t>||i on HI E0 (in dimension d — 3) 
is denoted in [10] with \\Dv ||, the symbol || || standing for the L 2 norm: see page 
40 of the cited reference. The equality ||i>||i = ||.Di>|| is essential for subsequent 
comparison between some estimates of ours and [10]). 

The exponential of the Laplacian. Let us put 

e tA v := e' m2 v k e k Vt G [0, +00), v G W ; (3.30) 
kez d 

the above series converges in the weak topology of W. 
From the Fourier representations, it is clear that 

e* A D' C W , e tA D' S0 c D' E0 Vt G [0, +00). (3.31) 

In the sequel, we consider any n G R. As well known, the map (t, v) 1— > e tA v sends 
continuously [0, +00) x H™ into H"; this map is a strongly continuous semigroup on 
the Hilbert space H™ , with generator A [ H™+ 2 : H[, t+2 -> H™. Furthermore 

||e tA w||n ^ e - *||u|| n Vt G [0, +00), v G (3.32) 

(as follows immediately from the Fourier representations (3.30) of e tA , combined 

with the obvious inequality e~*' fc ' ^ e~* for all k G Zf). 

Now, let us consider two Sobolev spaces H™, H™~^, with v G (0, +00); then 

e lA !)6iJ, lle^ulU < /j„(*) e_t |Mln-i/ V* G (0,+oo),u G (3.33) 



\2et/ .. 2 



V2et7 ' zT 2 ' (3-34) 

1 for t > - ; 
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(see Appendix A). The map (t,v) h- > e tA t> sends continuously (0, +00) x H™ " into 
EI?- 

All the previous statements about e tA and the Sobolev spaces, including the relations 
(3.32)-(3.34), hold as well if we replace systematically H™, H™+ 2 , M r { ^ u with H£ , 
H™+ 2 , H^T"; this is basically due to the relations AD' E0 = D' E0 and e tA D' E0 C D' E0 . 

NS equations. A setting for the incompressible NS equations on T d consists of 
three Banach spaces F + ,F,F_ of divergence free vector fields on T d , such that all 
conditions (P1-P6) of Section 2 are fulfilled by these spaces, by the operator 

A:= A \ F + : F + ^ F_ (3.35) 

and by the bilinear map 

9 : F x F -> F_ , (v, w) i-> 9(v, w) := -&(vmdw) . (3.36) 

In the above, v»dw is the vector field on T d with components (v»dw) r := Y^t=i v s d s w r ; 
the function spaces are chosen so that, for v,w e F, the pointwise products v s d s w r 
are well defined and y(f,w) belongs to F_. 

The initial value problem (2.12) corresponding to the above maps A, 9 and to any 
datum u G F takes the form 

Find u e C([0, T), F) sucii tiat 

u(t) = e tA u - f ds e (t - s)A £{u{s).du{s)) for all t G [0,T) , (3.37) 
Jo 

and is related, in the sense already indicated, to the Cauchy problem du/dt = Au — 
£(u»du), u(0) = Uq. One recognizes the NS evolution equation for an incompressible 
fluid (in units in which the density and the viscosity are 1, and assuming no external 
forces) . 

In [9] we considered for arbitrary d the setting F := H E0 , F± := M^f 1 (n G 
(d/2,+oo)), putting the emphasis on all the quantitative aspects and, in partic- 
ular, on the accurate estimation of the constants N, K in Section 2; this allowed to 
get fully quantitative estimates on the time of existence for the solution of (2.12), 
and gave a framework to evaluate a posteriori the distance between any approximate 
solution and the exact solution. 

As anticipated in the Introduction, here we will take a similar attitude for the setting 

F + :=H 2 r, F:=H E() , F_ := H E0 W , (3.38) 

mainly for d = 3 and 1/2 < uj < 1. Again, our evaluation of N, K will be rele- 
vant in connection with the time of existence for (3.37), with the estimates about 
approximate solutions, and so on. 
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4 An H 1 setting for the semigroup of the Lapla- 
cian 

For the moment we assume 

de {2, 3, ....}, we (0,1) (4.1) 

and take F_, F, F + as in (3.38), A as in (3.35); the inclusion (PI) F ^ F ^ F + 
is evident. 

4.1 Proposition, (i) A generates a strongly continuous semigroup on = F_; 
this is the restriction to [0, +oo) x of the map (t,v) ^ e tA v (see Eq. (3.30)). 
(ii) The map (t, v) \— > e tA v gives as well a strongly continuous semigroup on H^ = F; 
one has 

||e* A u||i <e-*||u||i Vt G [0,+oo),n G l£ . (4.2) 

(Hi) The map (t,v) i— > e tA n sends continuously (0, +oo) x = (0, +oo) x F_ into 
H^ = F; one /ias 

||e tA ^||i ^ ^ w (*)e _t ||u||- w Vt G (0, +oo), n G E£ W , (4.3) 

— e* /or 0<t( — — , 

\2et J 2 (4.4) 

(Note that ^(t) = 0{\/t^ 1+u) ^ 2 ) fort — > + ; £/iis power ofl/t is integrable due to 
the assumption < u < 1). With /i^ as above, one has 

sup / ds e~ s fi w (t — s) < +oo . (4.5) 

te[o,+oo) Jo 

(iv) Statements (i-iii) indicate that the set (F + , F, F_,A) fulfills all conditions (P2- 
P5) of Section 2. The function \i and the constants B,N G (0, +oo) mentioned 
therein can be taken as follows: 

/i — /lu as in (4-4), B — 1, N = N w := any majorant of the sup (4-5). (4.6) 

Proof, (i-iii) Essentially, one specializes to the present case all statements of Section 
3 (and Appendix A) about (e* A ) and Sobolev spaces of arbitrary order. In particular, 
the relations (4.2)-(4.4) follow from the relations (3.32)-(3.34) for the spaces H™ and 
H""", with n = 1 and v — 1 + uj. With this choice of u, the function fi u of Eq. 
(3.34) becomes the fj, u of Eq. (4.4); the finiteness of the sup (4.5) is easily inferred 
from the explicit expression of by an argument very similar to the one employed 
in [7] to prove Proposition 7.2, item (iv). 

(iv) Obvious. □ 
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5 An H 1 setting for the NS quadratic nonlinear- 
ity 

Throughout this section we assume 

de{2,3,...}, we^-Uoo). (5.1) 

Our aim is to show that the NS bilinear map (v,w) h- > < S > {v,vS) := — £(v»dw) is well 
defined from F x F = H* x H* to F = H~^, and to determine a constant K w 
such that || J ) (w, w)||-cj ^ 1 1 ^ 1 1 1 1 1 1 1 1 - 

We will reach this result through a number of intermediate steps, most of them 
relying on the function 

K.„: Z? - (0, +oo) , k~ Wk) := £ i _ ; (5.2) 

in the above one should intend 

Z^:=Z d \{0,A;}, (5.3) 

a notation to be employed systematically in the sequel. 
In Appendix C we show the following. 

(i) The sum defining K. w (k) (in principle, existing in (0, +oo]) is in fact finite, for 
each k G Zg. 

(ii) One has sup fe6Z d lC u (k) < +oo. 

(iii) Explicit lower and upper bounds can be given for JC W , to be used for the practical 
evaluation of this function and of its sup. 

The function K u first appears in the forthcoming Proposition, dealing with pointwise 
multiplication in the Sobolev spaces (of scalar functions) H", H^; this Proposition 
will be subsequently employed to estimate the NS bilinear map. 

5.1 Proposition. With d and uo as in (5.1), the following holds. 

(i) Let z G H£ , v G H], so that zv (the pointwise product of two L 2 functions) is 
well defined (and in L 1 ). This product has the additional features 

zv G L 2 , \\zv - (zv)\\ L 2 ^ T^llziyi^Hi (5.4) 

for a suitable constant K w G (0, +oo) 7 independent of z,v (recall that ( ) is the 
mean). 

(ii) The above inequality is fulfilled for all z,v by any constant such that 

'sup JC w (k) < Ku, (5.5) 



(2vr) d / 2 
with K.u as in Eq. (5.2). 
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Proof. Let us consider the Fourier coefficients (zv) k of zv, for k G Z d ; up to a factor 
l/(27r) d / 2 , these are obtained taking the convolution of the Fourier coefficients of z, v. 

^ V)k = J2^W E ZhVk ~ h = J2^W E ZhVk ~ h (5 ' 6) 
v ; hez d v ; hez$ k 

(the sum can be confined to Z^, due to the vanishing of z ,v ). From here we get 

< T^w^ E NK-^I ( 5 - 7 ) 

Now, Holder's inequality | J2h a h b h\ 2 < (Eft M 2 ) (Efe N 2 ) S ives 

IMfcl^^MW*) , (5.8) 

m*) := E ifei^ifc_ fe |2 as in ( 5 - 2 )> w := E i^ri^i 2 ^ - h\ 2 \vk-h\ 2 ■ 

Eq. (5.8) implies 



feez 



o 



^ f2^F ( SUP d MA;) ) E < ^ E ( 5 - 9 ) 



feeZ o 7 fcezg fcezjj 



(2tt) c 

where 1^ is any constant such that (^)~ d ^ 2 yJsup keZ d K u (k) ^ K^, as in (5.5). 
To continue, let us observe that the definition of V implies 

5X*o = E \hr\z h \ 2 \k - h\ 2 \ Vk _ h \ 2 

fcezg (fc,fc)ezgxzg,fc^/i 
whence, with a change of variable £ = k — h, 

y £v(k)= I^NWN 2 = INI«II«II?- (5-io) 

From Eqs.(5.9) (5.10), we get 
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The finiteness of ^2 keZ d \(zv )k\ 2 , and thus of J2kez d \( zv )k\ 2 , ensures zv G L 2 . Noting 
that zv — {zv) = X]feez ( t j(' 2 ' i; )fc e fc we can reexpress (5.11) as 

\\zv-(zv)\\ 2 L2 ^K 2 Jz\\l\\v\\l, (5.12) 

whence the thesis. □ 

5.2 Proposition. Let d, uj be as in (5.1), and denote with K u any constant 
fulfilling Eq. (5.5). Furthermore, let 

v G Ml , weM]; (5.13) 

then (i)(U) hold. 

(i) Consider the product v»dw (whose components {v»dw) r := J2^ =1 v s d s w r are well 
defined and belong to L 1 , being sum of products of the L? functions v s , d s w r ). One 
has 

v.dw G , \\vmdw\\- u s= iiT w ||u||i||«;||i . (5.14) 

(ii) Furthermore, one has 

2(v.dw) G H- w , \\£(v.dw)\\_u «C tfJuWHi . (5.15) 

(Due to the above results, condition (P6) of Section 2 is fulfilled with F := H^ , 
F : = I!-", 7 : F x F -> F , 9(v, w) := -&(vmdw) and K := Kj. 

Proof. Let us observe that (ii) follows immediately from (i) using (3.26) with 
n — > — uj and v — > v»dw. Hereafter, we give the proof of (i); after putting for brevity 

p := v»dw , (5.16) 

we proceed in several steps. 
Step 1 . One has 

p r G L] := L 1 fl D' (5.17) 

for all r G {1, d} (i.e., p G := L 1 fl W Q ). We have just noticed, as a comment 
in the statement of this Proposition, that p r G L 1 ; furthermore, p r has zero mean 
since 

/ p r dx = 5/ v s d s w r = - V / (d s v s )w r = (5.18) 

(here, we have used an integration by parts and the assumption divu = 0). 
Step 2. For all r G {1, d} one has 



P r eH;" , \\ P r \\^ ^ KM\ U 



\ 8=1 



15 



We know that p r G L\\ by the familiar duality between and H UJ , 



zp'etfVzeHZ, sup Lfr^jj <+00 . (5 20) 



if the above conditions are fulfilled, we further have 



(5.21) 



Keeping in mind the above statements, we consider any z G H" and the function 



zp 



= ^ z v s d s w r . 



(5.22) 



s=l 



For each s we have the following: 

(i) (9 s w r G L 2 , due to the assumption w G H„; 

(ii) G L 2 : this follows using Proposition 5.1, which also gives 

\\ZV S - (ZV S )\\ L 2 ^ ifa;||2;|| w |K||i . 

From zv s , d s w r G L 2 for all s, it follows 

zp r G L 1 ; 

now we estimate the integral of zp r , starting from the following remark: 



/ zp r dx = y { 



{zv s )d s w r 



(5.23) 



(5.24) 



(5.25) 



= V / (*u a - (^ s ))<9 s u? r d:r + YW) / 9 fl «; r da; = V / ( 



(zv s - (zv s ))d s w r 



(J T dd s w r dx = 0, since this is the integral of a derivative). From (5.25), from the 
Holder inequality and (5.23) we get 



zp r dx 



^5 ifJkllo: : 



^^2\\zv s - (zv s )\\ L 2\\d s w r \\ L 2 ^ \\d s w r \\ L 2 



s=l 



s=l 



.EMI, 



it) 



|2 

II 2 



^JklUMIi 



1 



£ \\e,w* HI, . (5.26) 
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Now, using (5.26) with (5.20) (5.21) we conclude that p r is actually in H^, and 
||p r ||-£j admits the bound (5.19). 
Step 3. One has 

peH ( ; u , IHU xJMHMIi (5.27) 

(so, (5.14) holds and the proof is concluded). From Step 2, we know that each 
component p r is in H~ u ; furthermore, the estimates (5.19) imply 

Ibll 2 -. = E Ibl 2 . < K\Hl E ll^||| 2 = A>||?H|? , (5-28) 

r=l r,s=l 

whence the thesis. □ 



6 Putting things together: an H setting for the 
NS initial value problem, in dimension d = 2, 3 

Let us recall that 

F + :=H£", F:=I4, F_ := ; A := A \ : M 2 ^ -> ; (6.1) 

3> : F x F -> F_ , (u, iu) i-> 3>(u, w) := -£(u.du;) . 

The above Sobolev spaces live, for the moment, on a torus T d of any dimension 
d G {2, 3, ..}; our previous results about this setting can be summarized as follows. 
Condition (PI) of Section 2 is fulfilled by the triple F + ,F,F_; according to Propo- 
sition 4.1, A fulfills conditions (P2-P5) if uj G (0, 1). On the other hand, accord- 
ing to Proposition 5.2, is well defined and fulfills condition (P6) of Section 2 if 
uj G {d/2 — 1, +oo); so, the requirements of both Propositions 4.1, 5.2 hold simulta- 
neously if 

^ - 1< u < 1 . (6.2) 

Of course, (6.2) can be fulfilled for some uj only if d/2 — 1 < 1, i.e., for d = 2 or 
d — 3; so, (6.2) holds if either 

d = 2, we (0,1) (6.3) 

or 

d = 3, we(j,l). (6.4) 

Summing up (and recalling the statements of Propositions 4.1, 5.2 on the function 
ji and the constants B, N, K), we have the following result. 
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6.1 Proposition. Let d,cu be as in Eqs. (6.3) (6.4)- Then, the set (F + , F, F_, A, 1 ?) 
defined by Eq. (6.1) fulfills all conditions (P1-P6) of Section 2. The function /i and 
the constants B, N, K mentioned in Section 2 can be taken as follows: 

\i = [1^ as in (4-4), B — 1, N = as in (4-6), K = as in (5.5) . (6.5) 

The previous proposition allows to apply the full machinery of [9] to the NS initial 
value problem (3.37). Hereafter we consider, in particular, the condition for global 
existence presented in [9] and summarized in Proposition 2.1 of the present work. 

Global existence for NS with small initial data, when d = 3. Let us keep 
the definitions (6.1). Global existence for the NS initial value problem (3.37) is well 
known for any datum in F = H* , if d = 2; so, we pass to the case d = 3. Let us recall 
that, in the three dimensional case, we have the equality (3.29) ||i>||i = 1 1 curl v 1 1 x,a for 
all v G M^ . On the grounds of Propositions 2.1 and 6.1, we can state the following. 

6.2 Proposition. Let d = 3, lu £ (1/2,1) as in (6.4), and F, etc., as in (6.1). 
The solution u of the initial value problem (3.37) is global if the initial datum u £ 
F = H^ is such that 

IMI. < ^ I (6-6) 

furthermore, the solution fulfills the bound ( of the type (2.14)) 

||u(t)||i ^ A^KJVJKIIOe-'Klli f or t G [0, +00) , (6.7) 

where X G C([0, 1], [1, 2]) is the increasing function defined by (2.15). More roughly, 
we have 

||u(i)||i < 2e"* ||u ||i for t G [0, +00) . (6.8) 



To get a fully quantitative estimate, let us put 



w := ^ (6.9) 

(see the forthcoming Remark 6.4 (iii) about this choice). Then, computing numeri- 
cally the function of t indicated below, we see that 

sup / dse~ s ix 7/w (t - s) < 1.70 := N 7/w (6.10) 
te[o,+oo) Jo 

(in fact the above sup is a maximum, attained at point t G (0.21, 0.22); for these and 
other numerical computations one can use, e.g., the MATHEMATICA package). 
Furthermore, one has (see Appendix D) 

27.94 < sup lC 7/w (k) < 32.23 (6.11) 
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(here we are mainly interested in the upper bound 32.23; the lower bound 27.94 is 
reported only for an appreciation of our uncertainty about sup/Cy/io). This implies 



0.335 < j—^ /sup /C 7/10 (£;) < 0.361 := K 7/w . (6.12) 



The values of N 7 / 10 and K 7 / 10 imply 

1 



4N 7 / W K 7 / 10 



> 0.407 ; (6.13) 



In conclusion, Proposition 6.2 with the choice u = 7/10 and the previous evaluations 
of the related constants yield the following. 

6.3 Corollary. With d = 3, the solution u of the initial value problem (3.37) is 
global if the initial datum uq G F = H* is such that 

\\u \\i ^ 0.407 ; (6.14) 

furthermore, for all t e [0, +oo) ; the solution fulfills the bounds 

IKOIIi < e ~* l|Mo111 ^ 2e ~* l|Mo111 ' (6 - 15) 

with X as in (2.15). 

6.4 Remarks, (i) According to the remark following Proposition 2.1, the global 
solution u can be constructed by a Picard iteration. A fine analysis of the iteration, 
based on the the smoothing properties of (e* A ), shows that the function (t,x) h- > 
u(t)(x) = u(x,t) is in fact C°° on (0, +oo) x T 3 , where it satisfies the NS equations 
in the classical sense; on this point see, e.g., the proof of Theorem 15.2(A) in [4]. 
Taking into account these facts and the exponential time decay of ||-u(t)||i, we see 
that u is in L°°([0, +oo), L| ) n L 2 ([0, +oo),H^ ) and fulfills the NS equations in the 
weak sense; so, u is a weak solution of the NS equations in the sense of [11]. Again 
by the exponential time decay, it is u e £ 4 ([0, +oo), H^ ); this suffices to infer that 
u is a strong NS solution in the sense of [11] (see the Remark 3.3 on pages 22-23 of 
this reference). 

(ii) We already mentioned that our criterion (6.14) improves the condition for the 
existence of a global (strong) solution recently proposed in [10]; in our notations, 
the condition of the cited paper reads 

||u ||i ^ 0.00724 (6.16) 

(see page 43 of [10]; 0.00724 is the approximation with three meaningful digits of 
the quantity indicated therein by Ry)- 
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(iii) Of course, in place of to = 7/10 one could consider for cu other choices in the 
interval (1/2,1), each one yielding a bound of the type (6.6) for global existence; 
the best estimate of this type would be obtained maximizing 1/(AN U K U ), for uo in 
(1/2,1). However, a few experiments we did with uo ^ 7/10 seem to exclude a 
significant improvement of the bound (6.14). 

7 Some possible developments. 

In [9], we presented the following abstract result: 
7.1 Proposition. Assume the following: 

(i) (F + , F, F_,A,T) is a set with the properties (P1-P6) of Section 2; fi,B,N,K 
are the function and the constants mentioned therein. 

(ii) u G F is a datum for the initial value problem (2.12). 

(iii) u ap is an approximate solution of problem (2.12) with domain [0, T) (T G 
(0, +oo]j, a growth estimator!) and an error estimator £ : these expressions indicate 
three functions u ap G C([0, T),F) andV,£ G C([0,T), [0, +oo)) such that 

\\u ap (t)\\^V(t) , (7.1) 

\\u ap (t)-e tA u - [ dse^ A 7(u ap (s),u ap (s))\\^£(t) forte[0,T). (7.2) 
Jo 

(iv) There is a function R G C([0,T), [0, +oo)) fulfilling the "control inequality" 

S(t)+K ds n{t - s)e- B{t - s) (2V(s)R(s) + R 2 (s)^j <: R(t) for t G [0, T) . (7.3) 

Then, problem (2.12) has an (exact) solution u G C([0,T), F) , and 

\\u(t) - u ap (t) || ^ R(t) for t G [0, T) . (7.4) 

The previous proposition allows to make predictions on the interval of existence of 
the solution u of (2.12), and on its distance from u ap , using some information (T> 
and £) pertaining to u ap only; in this sense, we have an estimate for u from an a 
posteriori analysis of u ap . 

In [9], we presented some applications of this result to the NS equations on T d , 
taking for F ± , F some Sobolev spaces of sufficiently high order and considering, for 
example, the Galerkin approximate solutions. 

The analysis performed in the present work allows to apply Proposition 7.1 to the 
NS equations with F + := H|"^,F := H^ ,F_ := E£ W and (d = 2, oj G (0,1) or) 
d = 3, u G (1/2,1). 
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Our estimates on some related functions and constants (e.g., \i = [i w and K = K u 
considered in Sections 4-6) allow a fully quantitative implementation of the method. 
It should be pointed out that the availability of accurate information on /j,, B, K, etc., 
is essential for an efficient application of the control inequality (7.3); for example, 
the time interval [0, T) on which a solution R of (7.3) exists depends sensibly on 
these data. We plan to illustrate elsewhere some applications of Proposition 7. 1 to 
the NS equations on T 3 , in the framework of the above spaces F±, F. 

A Appendix. On the semigroup of the Laplacian 

After recalling the definition (3.30) 

e tA v := e ~ m2v kZk V* G [0, +oo), v G W , 
fcez d 

let us fix n G R, v G (0, +oo). Hereafter, we prove the following result. 

A.l Proposition. Let t G (0, +oo) ; v G H™ _1/ . Then, we have the relations 
(3.33), (3.34) 

e tA v G H™ , ||e tA t;|| n < /i^K* IMU-, , 



K{t):={ ^et) ^ 2 



1 fort> - . 

J 2 

Proof. First of all, since v G W we have e tA v G D' . To go on, let us observe that 

\k\ 2n \(e tA v) k \ 2 = \k\ 2u e- 2t W 2 \k\ 2n - 2 »\v k \ 2 (A.l) 

< ( sup Ik^e'^ 2 ) f \k\ 2n - 2v \v k \ 2 ) < ( sup U^oj) \\v\\ 2 n _ u , 

Vjp) := re~ 2ti > . 
An elementary computation gives 

f uJ^) = (^Y for0<t^, 
sup U v 0) = }2t) Wt) v 2 (A.2) 

*e[i,+oo) y U vt {l) = e~ 2t for t > - . 

Returning to (A.l), we infer that e tA v G HI™, and 

|e* v\\ n ^ a j sup U ut (d) \\v\\ 

n—v i 

(A.3) 

J?6[l,+oo) 
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expressing the above sup via Eq. (A. 2) and isolating a factor e * we easily get Eqs. 
(4.3) (4.4). □ 



B Appendix. Some results on the convolution. 

The results we present here are used in the next Appendix to establish some facts 
about the function K. w , defined on Z d via Eq. (5.2); the cited equation tells us that 
JC w (k) is a convolutionary sum (apart from a technical detail, i.e., the elimination 
of and k from the summation domain Z d ; we return on this later on). 
In this Appendix we give some general results on convolutions, which are stated 
independently of the subsequent applications to /C w . Of course, the convolution of 
two functions /, g : Z d — > [0, +oo) is defined by 

/ * g : Z d - [0, +oo] , fc -> (/ * g)(k) := £ f(h)g(k - h) . (B.l) 

hez d 

We first consider the case d — 1; our main statement for this case is contained in 
Proposition B.2, which is preceded by the following definition. 

B.l Definition. Consider a function 

/:Z-[0,+oo], k^f(k). (B.2) 

/ is even if 

f(-k) = f(k) for allkeZ; (B.3) 

/ is unimodal if it is nondecreasing on — N = {0, — 1, —2, ...}, and nonincreasing on 
N = {0,1,2,...}: 

f(£) ^ f{m) for £,m G — N and £ ^ m , 

f(£) ^ f(m) for£,meNand£^m. (B.4) 

Obviously enough, the unimodality of an even / is equivalent to any one of these 
conditions: / is nonincreasing on N, or 

f(£) ^ f(m) for £,meZ and \£\ < \m\ . (B.5) 

B.2 Proposition. Let us consider two even, unimodal functions 

p,g:Z^[0,+oo) , (B.6) 

and put 

s :=p*q : Z -> [0, +oo] , k i-> s(k) = ^p{h)q(k - h) . (B.7) 
Then, s is itself even and unimodal. 
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Proof. (Adapted from the proof of Proposition 4.5.5 in [5]; the cited result refers 
to continuous convolutions, with integrals over R in place of sums over Z). 
To prove that s is even we write, for any k e Z: 

s (~ k ) = $>(%(-^-^) = $>(-J'M-^ + j) = $>(J'M^-J) = s(k) (B.8) 
hez jez jez 

(the second equality relies on the change of variable h — —j, the third one holds 
because p, q are even) . 

Proving the unimodality of s is less trivial. It suffices to prove that s is decreasing 
on N, which can be written as follows: 

s{£)^s{£+l) for^GN; (B.9) 

in the sequel, we fix any £ and derive the thesis (B.9). First of all, we write 

— 1 oo 

s (£) = $>(%(£ -h)=Y, ~ h ) + Y,p^ £ - h ) ; 

h£Z h=-oo h=0 

now, putting h = —k — 1 in the first sum, h = k in the second one and using 
p(-k - 1) = p(k + 1), we get 

oo 

«W = ( P( k + ^ + ^ + !) + p(*0?(* - fc) )■ (B.10) 
Similarly, writing 

oo 

s(£+l) = J2p(h)q(£-h + l)= p(h)q(£-h + l) + J2p(hW-h + l) , 

h&Z h=-oo h=l 

putting h = —k in the first sum, h = k+1 in the second one and using p(—k) = p(k), 
we get 

oo 

s(£+l) = J2( P(k)q(£ + k + 1) + p(k + l)q(£ - k) ). (B.ll) 

fc=0 

Now, consider any k e N. Then = fc ^ + 1 = \k + 1| and \£ — k\ ^ \£\ + 
= £ + A;^£ + A; + l = |£ + A; + l|so that, by Eq. (B.5) with / = q and / = p, 

p(k)^p(k + l), q{£-k) ^q(£ + k + l) . (B.12) 

Thus ^ \p(k) - p(k + l)][q(£ - k) - q(£ + k + 1)] = p(%(£ - k) + p(ib + 1) 
q(£ + k + 1) - p(%(£ + fc + 1) - + l)g(£ - k), i.e., 

p(Jfe + l)q(£ + k + 1) + p(%(£ - fc) ^ p(k)q(£ + k + 1) + p(k + l)g(£ - fc) . (B.13) 
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Now, from (B.10) (B.ll) (B.13) we get the thesis (B.9). □ 

Let us extend the previous considerations from the one-dimensional to the d-dimensional 
case, for arbitrary d. 

B.3 Definition. Consider a function 

/:Z d -[0,+oo] , &-/(*;); (B.14) 
for each r G {1, d} and k' = (k±, fc r _i, k r+ i, ka) G Z d_1 7 pwt 

/ rfc # : Z -> [0, +00] , k r i-> f rk '(k r ) ■= f(h, fc r _i, fc r , fc r+ i, fc d ) . (B.15) 

VFe call f even (resp., unimodal) in each variable if, for any r G {l,...,d} and 
k' G Z d ~ x , the function f r y is even (resp., unimodal) in the sense of Definition B.l. 

We note that / is even in each variable if and only if 

f(Rrk) = f(k) for each r G {1, d} and k G Z d , (B.16) 

R r k := [ki, k r _i, —k r , k r +±, k^) . 

If / is even in each variable, recalling Eq. (B.5) we see that the unimodality of / in 
each variable is equivalent to 

f(f) > f(m) if £,me Z d and < \ mi \, \£ d \ ^ \m d \ . (B.17) 

B.4 Proposition. Let us consider two functions 

p, q : Z d -> [0, +00) (B.18) 

which are even and unimodal in each variable, and put 

s :=p*q : Z d -> [0, +00] , k ^ s(k) = ^ p(h)q(k - h) . (B.19) 

TTien, s is rise// even and unimodal in each variable. 

Proof. Let us fix r G {1, d}, k' = (k± : fc r _i, fc r +i, fed) G Z d_1 , and consider 
the function s r /j : Z — > [0, +00], defined following Eq. (B.15); we must prove that 
s r y is even and unimodal. 

To this purpose we note that the definition s := p * q implies 

s rk >(k r ) = ^2 (Prh> * q r ,k'-h')(K) for each k r G Z , (B.20) 



24 



where the functions p r h',Qr,k'-h' '■ Z — > [0, +00) are defined as well following Eq. 
(B.15). Each of the functions p r h>, q r ,k'-h' is even and unimodal; so, due to Proposi- 
tion B.2, their convolutions p r h> * q r ,k>~h> are even and unimodal. A sum of even and 
unimodal functions has the same properties, so our thesis about s r k' is proved. □ 

Our last statement about convolutions is obvious, and mentioned only for subsequent 
citation. 

B.5 Definition. A function 

/:Z d ^[0,+oo] , fc-/(fc) (B.21) 

is symmetric if, for each a in &d (the group of permutations of {1, ...,d}) and each 
k = (ki, k d ) G Z d , it is 

f(P a k) = f(k) , P a (k) := (k a{1) , k a{d) ) . (B.22) 

B.6 Proposition. Ifp,q : Z d — > [0, +00) are symmetric functions, their convolu- 
tion p * q is itself symmetric. 



C Appendix. The function K 



Let d G {2, 3, ...}. In the present Appendix we prove a number of properties of the 
function /Q, on Z d , defined by Eq. (5.2); some of these properties were mentioned 
in Section 5. Some of our results rely on the following two Lemmas. 

C.l Lemma. Let us consider two radii p, p x such that 2\fd ^ p < p\ ^ +oo ; and 
a nonincreasing function x G C([p — 2\fd, p x ), [0, +00) ). Then, 

E x(\ h \) ^ WJ7^\ / dt(t + Vd) d ~\(t) (C.l) 
(the two sides of the inequality being, possibly, +00). 

Proof. See the Appendix C of [9]. □ 

C.2 Lemma. Fix v G (d, +00) and consider, for any A G (0, +oo) ; the sum 

AS„(A):= E ]4; ( C -2) 
hez d ,\h\^x+2Vd 

this is finite, and admits the bound 

O n d/2 *zl / A _ 1 \ d d/2-l/2-i/2 

*MD < «.(A) == ^ £ ( , ) (y ! i _ 1)A „-,- 1 • (C3) 
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Proof. We apply the previous Lemma, noting that 

AS„(A)= X(H), p:=A + 2v/rf, Pl :+oo, X (t) := i . (C.4) 

ftez d ,p^|fe|<pi 

The inequality (C.l) gives 

^ (A)< fP)A dt — F — ; (C - 5) 

now, writing the expansion (t + v^)** -1 = X^=o ( ^ j ^ ) t l d d l 2 ~ 1 l 2 ~ 1 / 2 and inte- 
grating term by term, we see that 

r.h.s. of (C.5) = 8$ U (X) . (C.6) 

The relations (C.5) and (C.6) give the thesis. □ 
From now on, we make the assumption (5.1) 

u G - l,+oo) ; 
as in Eq. (5.2), we consider the function 

K u :keZ d » }C u (k) := 



\h\ 2 "\k-h\< 



-■ok 



As noted in the comment following (5.2), the sum defining /C w certainly exists in 
(0, +00] ; hereafter we show its finiteness, with many other properties of the function 
under investigation. 

C.3 Proposition, (i) For each k G Z d , it is K u {k) < +00; furthermore, for any 
"cutoff" A G {\k\, +oo) ; one has 

X u (k, A) < JC w (k) < 3C W (A;, A) + <S0C w (A;, A) (C.7) 
*UM):= £ ' (C - 8) 

feeZ ( '5 fc ,|fe|<A+2v / d 

"^^iwM • J(2 U + i-o(a-i*|)^- (G9) 



i=0 



T/ie function k \— > lC w (k) is even in each variable and symmetric, i.e., 

K, u (Rrk) — K. w (k) for all r G {1, ... ,d} and k e Z d , (CIO) 
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IC^k) = K u {k) for alia e& d and k eZ d (C.ll) 

(where, as in Appendix B: R r k := (k±, — k r , ka), &d (ire the permutations of 
{1, d}, and P a k := (k a(1) , k a{d) )). 
Furthermore, one has 

M*) + ^ + ^ < M^) + + (C.12) 

fMezi «nd | fci| > ^ |^| . 

(Hi) It is sup fc6Z <* K, w {k) < +00. For any a G {1, 2, ...} ; we have the bound 

sup /C w (&;) < max(5 u (o),7J il) (a)) , (C.13) 



^(a) := max /C w (Jfe) , /(a) := {A; G | ^ h ^ fc 2 ... ^ k d ^ a} , (C.14) 

fceJ(a) 

T^a) := /CCO, 0, a + 1) + + . (C.15) 

Proof, (i) Of course, the bounds (C.7-C.9) to be proved imply finiteness of the sum 
/C w (fc); let us derive these bounds, for any fixed fc G Z[[ and cutoff A G (|fc|,+oo). 
First of all, we write 

% u {k, A) < K u (k) = XUk A) + AX^k, A) , (C.16) 

where %u>(k, A) is the finite sum in (C.8), and 

*MM)== £ jSjSjjbiF <C - 17) 

/iGZ d ,|/i|^A+2v / rf 

(note that ^ A + 2^ implies h 7^ 0,/c). In principle, AX u) (k, A) G (0, +00]; 
hereafter we will prove that 

AX u (k,X) < <WC W (M) (C.18) 

where 5X Ld (k, A) G (0, +00) is defined by (C.9); this will give the inequality (C.7) 
that, with the finiteness of 5% w {k, A), also implies lC u (k) < +00. In order to prove 
(C.18), we put 

p:=^^, q:=u + l; (C.19) 
ui 

then p, q G (1, +00) and 1/p + 1/q = 1 so that, by the Holder inequality, 
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wx( E E 

heZ*,\h\^\+2y/d h£Z d ,\h\^\+2Vd 

= ( E JhP^) ( ^ |jfe -^+2) • (C - 20) 

heZ d ,\h\^\+2Vd hdZ d ,\h\>\+2Vd 

Let us consider the two sums appearing in the last passage of (C.20). A change of 
variable h' — k — h in the second one gives 

Ifc _ /J 2^+2 ~~ E/ |/ 1 /|2oj+2 ' (^-21) 

/ieZ d ,|/i|5sA+2 v / d /i'eZ d ,|/ l '-fc|^A+2v / d 

On the other hand, the inequality \h' — k\ ^ A + 2\fd implies \h'\ = \{h' — k) + k\ 
^ \ \h' — k\ — \k\ \ ^ X + 2\fd — \k\; so, the domain of the last sum is contained in the 
domain {\h'\ ^ A + 2\fd — and we conclude 

E/ Ifc _ /i|2w+2 ^ |/ 1 /|2oj+2 ' (C.22) 

/ieZ d ,|/i|>A+2v / d fe'eZ d , \h'Xz\+2\fd-\k\ 

Concerning the first sum in (C.20), it is obvious that 

E \h\ 2 ^+ 2 ^ \h\ 2 "+ 2 (C.23) 

h£Z d ,\h\^\+2Vd hez d ,\h\>\+2Vd-\k\ 

(since the right hand side is a sum on a larger domain). We return to (C.20), and 
insert therein the bounds (C.22) (with h! renamed h) and (C.23); the conclusion is 

/ x-^ 1 \l/P+l/9 

A^(Jb,A)<( £ p^) 

feeZ d ,|fe|^A+2v / d-|fe| 

E ^4^ = AS 2 . +2 (A-|A;|)^5S 2w+2 (A-|A;|), (C.24) 

/iez d ,|fe|^A+ 2v / d-|fc| 

where the last two relations follow, respectively, from the definition (C.2) of AS^A) 
and from the bound (C.3), here applied with v — > 2u> + 2 and A — > A — (note 
that 2uj + 2 > d, by the assumption (5.1)). On the other hand, explicitating the 
definition (C.3) of 5§ v we see that 

5S 2w+2 (A - \k\) = 5X„(k, A) as in (C.9); (C.25) 

with (C.24), this yields the thesis (C.18). 
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(ii) We want to show the properties (C.10-C.12) of KL W . To this purpose, let us define 

p : Z d ^ [0, +00) , k h-> p(k) := { * * f °| (C.26) 

9 : [0, +00) , k h-> «,(*) := | j^ 1 ' * * f J (C.27) 
Then, for fceZf, 

(p * ?) (*0 = - /i) = J] p(/i)g(A; - h) + p(0)q(k) + p(fc)g(0) ; 

hez* he zg fc 

explicitating p and q, we can rephrase this as 

(p* q )(k)=JC u (k) + -^ + ^ ; for fceZ„ . (C.28) 

The functions p, g are even and unimodal in each variable, as well as symmetric (in 
the sense of Appendix B); by Propositions B.4, B.6, the same properties hold for 
their convolution p * q : Z d — > [0, +00). So, we have 

(p*q)(R r k) = (p*q)(k) for all r G {1, d} and k G Z d , (C.29) 

(p * g) (i^Jfe) = (p*q) (k) for all a G & d and fc G Z d , (C.30) 

(p*q)(k) s$ (p*q)(i) ifk,£eZ d and \k r \ ^ \£ r \ forr = l,...,d (C.31) 

(the last relation is the inequality (B.17) expressing the unimodality of / := p * q, 
with the replacement m — > k). For k G ZjJ, explicitating (p * q)(k) via Eq. (C.28) 
(and using the obvious relations \R r k\ = \P a k\ = \k\) Eqs. (C.29) (C.30) (C.31) 
yield, respectively, the conclusions (CIO) (C.ll) (C.12). 

(iii) Let us prove Eq. (C.13) (which, of course, implies the finiteness of sup fcgZ d K, w {k)). 
To this purpose, we choose any a G {1, 2, ...} and write 

Z d = J (a) U L(a) , (C.32) 
J(a) := {k G Z d \ \k r \ < a for all r G {1, 4} , 
L(a) := {k G Z d \ \k s \ ^ a + 1 for some s G {1, d}} . 

Hereafter we will derive the bounds 

M fc ) < S M for all k G J(a) , (C.33) 

K, u (k) < K,(a) for all fc G L(a) , (C.34) 
yielding the thesis (C.13). 
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Let k G J (a); applying a reflection to each negative component of k (if any), and 
then performing a suitable permutation, we can transform k into an element of the 
set 1(a) in Eq. (C.14); more formally, there is a map C = P a R rm ...R ri of into 
itself such that Ck G 1(a). Due to the results of (ii), K.^ is invariant under C . So, 

K u (k)=lC u (Ck)^S u (a) , 

and (C.33) is proved. 

Now, consider any k G L(a). Then |/c s | ^ a + 1 for some s G which 
obviously implies 

\k r \ > |4| for r = l,...,d , £ := (0, a_+_l, 0) ; 

place r 



from here and (C.12) we infer 



" ■ |/|2 ' /12a; 

2 



On the other hand, the definition of £ and the symmetry of /Q, give fC^(£) +l/\£\ 
+ l/\£\ 2 " = K u (0, 0, a + 1) +l/(a + l) 2 + l/(a + l) 2w = 7^ (a), so 

this result is even stronger than the desired relation (C.34). □ 



D Appendix. Evaluation of sup fcGZ 3 K>7/io(k). 

We specialize the results of the previous Appendix to the case 

d = 3, ^ = 75; (D-l) 
our aim is to justify the statement (6.11) 

27.94 < sup K,y w (h) < 32.23 . 

First of all, let us we write down for the sup of )C 7 /i the bound (C.13), with a — 1. In 
this case 1(a) = 1(1) = {(0, 0, 1), (0, 1, 1), (1, 1, 1)}, so <S 7/10 (1) = max (/C 7/10 (0, 0, 1), 

/C 7/10 (0,1,1), /C 7/10 (l,l,l)); furthermore U 7/10 (l) = /C 7/10 (0, 0, 2) +1/4 + 1/2 7 / 5 . 

Eq. (C.13) states that sup feeZ 3 /C 7 /i ^ max ^<S 7 / 10 (1), 7?. 7 /i (l)^ , i.e., 
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sup JC 7/10 (k) (D.2) 

< max (/C 7/10 (0, 0, 1), /C 7/10 (0, 1, 1), /C 7/10 (l, 1, 1), /C 7/10 (0, 0, 2) + 1 + J^) . 

To evaluate /C7/10 at the points indicated above we use Eqs. (C.7-C.9), with a cutoff 
A = 150. The results are: 

27.94 < /C 7/ io(0, 0, 1) < 32.23; 27.48 < /C 7/10 (0, 1, 1) < 31.77; 26.49 < /C 7/ i (l, 1, 1) < 30.78; 

25.69 < /C 7/10 (0, 0, 2) < 29.98; /C 7/10 (0, 0, 2) + J + ^ < 30.61 . (D.3) 

(Note that the reminder term £JC 7/ i (£;, A) of Eqs. (C.7-C.9) is 0(1/A 2 / 5 ) for A -> 
+00; such a slow decrease at infinity explains why the lower and upper bounds in 
(D.3) are not very close, even with the fairly large chosen cutoff A = 150.) 
From (D.2) and the upper bounds in (D.3) we conclude 

sup K 7/W (k) < 32.23 , 
fcezg 

as stated in (6.11). Of course, /C 7 / 10 (0, 0, 1) ^ sup fceZ 3 JC 7 / w (k), and the lower bound 
for the former in (D.3) implies 

27.94 < sup £7/10 (A;) ; 
fcezg 

our justification of (6.11) is concluded. 
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